Droplet detachment and bead formation in visco-elastic fluids 
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The presence of a very small amount of high molecular weight polymer significantly delays the 
pinch-off singularity of a drop of water falling from a faucet, and leads to the formation of a long-lived 
cylindrical filament. In this paper we present experiments, numerical simulations, and theory which 
examines the pinch-off process in the presence of polymers. The numerical simulations are found to 
be in excellent agreement with experiment. As a test case, we establish the conditions under which 
a small bead remains on the filament; we find that this is due to the asymmetry induced by the 
self-similar pinch-off of the droplet. 



The pinch-off of liquid drops has attracted consider- 
able attention in recent years 0, not the least owing to 
its enormous technological applications in biotechnology, 
micro-scale manufacturing, and spray technology 
More recently, the breakup of droplets of visco-elastic liq- 
uids has been studied both experimentally and theoreti- 
cally. The interest of considering visco-elasticity is both 
fundamental and applied. From the fundamental side, 
the addition of minute amounts of polymers has been 
shown to inhibit the finite-time singularity that happens 
at breakup. From the applied side, in many applications, 
such as fire fighting, inkjet printing, or pesticide deposi- 
tion on plant leaves, complex fluids have been used to 
control or modify drop sizes. 

For drop detachment in Newtonian fluids, the dynam- 
ics is governed by the mathematical properties of the sim- 
ilarity solutions which describe the thinning of the fluid 
neck. Most importantly for applications, these similar- 
ity solutions are always independent of initial conditions, 
and hence it is very difficult to control drop breakup by 
an external manipulation of parameters. Therefore, the 
opportunities offered by complex fluids are tremendously 
important. 

In this Letter, we consider the detachment process in 
the simplest possible visco-elastic case, that of a dilute 
aqueous solution of flexible polymers. The detachment is 
studied experimentally for a well-characterized polymer. 
This allows us to use a model for the elasticity due to 
the introduction of the polymers, and perform numerical 
simulations of the detachment process. The good agree- 
ment between experiment and numerics permits to study 
the detachment process in great detail. As a test case, 
we investigate the condition under which the so-called 
"beads on a string" structure forms [5| , shown in FigHJi- 
This intriguing structure is ubiquitous in complex fluids, 
having been observed for solutions of both flexible and 




FIG. 1: The droplet detachment process of an aqueous 
lOOppm PEO solution at different time steps t c — t relative to 
the time t c at which the filament is formed. The left half of 
the pictures are plots from numerical simulations, the right 
half are experimental photographs : a): nozzle radius R = 
1.5 mm (pictures are 1.2 x 9 mm) t c — t= 6, 2, 0, -3, -5 ms; 
b: R = 0.4 mm (pictures are 0.9 x 0.9mm) t c — t= 1, ms. 
The model parameters used for the numerical simulation were 
rj p = 3.7 x 10 _4 Pas for the polymeric contribution to the vis- 
cosity, a polymer timescale of A = 1.2 x 10~ 2 , and an elasticity 
parameter of 6 = 2.5 x 10 4 . The solvent viscosity (water) is 
n s — 1 x 10 _3 Pas and the surface tension is 7 = 6 x 10 -2 N/m. 



rigid polymers in both low-and high viscosity solvents, 
and various semi-dilute and concentrated wormlike mi- 
cellar surfactant solutions. However, depending on the 
experimental conditions, beads may or may not form. 

For the experiments, we focus on the particular case 
of a drop falling from an orifice (cf. FigQ, although the 
theoretical arguments we develop here are more general, 
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and also apply for instance to filament stretching or fil- 
ament thinning rheometry. Experiments were performed 
with aqueous solutions of polyethylenoxide (POE) with 
a molecular weight of 4 x 1 6 amu, in concentrations of 5 
to 2000 ppm. Nozzle radii R ranged from R = 0.25mm 
to R = 5mm and the droplets were generated using a sy- 
ringe pump in a quasi static mode. Pictures were taken 
with a high speed camera (1000 frames/s). 

We find that the initial stages of the thinning of the 
liquid column are well described by Rayleigh's linear the- 
ory [l| , which predicts an exponential growth of the most 
unstable mode. Thus a growing trough forms on the liq- 
uid bridge that separates the falling drop from the tap, 
whose shape is fairly symmetric around the minimum 
bridge thickness, as seen in the first panel of FigQJi and 
b. In the second panel of Fig^, the neck shape has al- 
ready turned asymmetric, while for the smaller pipette 
radius the shape turns directly to a uniform cylinder. 
For the asymmetric case, the remaining frames of Fig^, 
show how the asymmetric neck structure turns into a 
secondary bead. 

In order to model the detachment, we use the FENE- 
P (Finite Extensibility Non- linear Elasticity) model [6j, 
a simple polymer model which treats the polymers as 
Hookean springs of finite extensibility. The springs can 
be stretched by the flow, leading to elastic stresses due to 
the entropic elasticity of the polymer chains. To further 
simplify the description, we use a Id long-wavelength the- 
ory for the fluid motion 0, leading to a coupled set of 
equations for the local radius h(z,t) of the fluid column, 
the mean fluid velocity v(z,t), and the radial and axial 
components of the polymeric stress cr r (z, t), <r z (z, t). For 
the present problem, only the component in the axial di- 
rection <r z (z, t) is of interest, because this is the principal 
direction of stretching. 

The model parameters used in the simulation were 
taken from the rheological measurements of 0- To ob- 
tain an optimal description, we used a polymer relaxation 
time A = 1.2 x 10 _2 s, about 4 times as long as suggested 
by , from steady-state rheological measurements. Our 
comparison theory/experiment can be seen as a method 
of determining the timescale relevant for polymer-flow in- 
teractions in a highly non-equilibrium situation. Given 
the considerable simplification inherent in a single-time- 
scale approach, we believe the correspondence in time 
scales between the rheological method and the value used 
here to be quite reasonable. 

In Fig. ^ experiments are compared to numerical sim- 
ulations, finding excellent agreement for the two different 
pipette radii was obtained by adjusting a single parame- 
ter, the polymer relaxation time A. For the larger pipette 
radius of 1.5mm a secondary bead is formed (cf. FigQ] 
a)), for the smaller radius of 0.4mm the elastic thread is 
completely uniform. 

Fig. |2K shows a phase diagram for the existence of 
a bead. We find that for low polymer concentrations 
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FIG. 2: a)The phase boundary for bead formation, obtained 
by analyzing experimental runs at polymer concentrations 
of 5, 10, 20, 50, 100, 200, 500, 1000 and 2000 ppm and nozzle 
radii R = 0.25, 0.75, 1, 1.5, 2, 3, 4, 5mm (triangles). Between 
the dashed and the full lines a bead is formed, but subse- 
quently disappears, b): The asymmetry parameter a for the 
the 100 ppm runs measured at t c — t — 1ms. The vertical 
lines correspond to the transition lines in a). 



and large nozzle radius a bead is observed, while in the 
opposite corner, below the full line, no bead is formed. 
In the narrow strip between the dashed and the full lines 
a bead is formed, but is so small that it gets stretched 
out by elastic stresses and eventually disappears. 

From the data shown in Fig. 1 , it appears that the tran- 
sition to the asymmetric neck plays an important role for 
the formation of the secondary bead structure. We quan- 
tify the asymmetry using the parameter a 0: 



J k {Zmin ~\~ % ) tl (Zrnin % ) dz 



h(z r , 



(i) 



whose value is zero for a symmetric and one for a stair- 
case profile. The boundaries of integration were chosen 
at h (±z') = 2h min . In Fig|2> we show a cut through the 
phase diagram at constant polymer concentration. In- 
deed, we find that the asymmetry parameter a just before 
onset of the filament undergoes a rather abrupt transi- 
tion from zero to a nonzero value, exactly at the point 
corresponding to the boundary between the bead and no- 
bead regions of Fig|2]} in the phase diagram. Similar cuts 
at constant nozzle radius show the same behavior. 
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FIG. 3: The minimum neck radius versus time for R = 0.4mm 
(circles), R — 1mm (triangles), and R = 4mm (squares), at 
a concentration of lOOppm, where to is the extrapolated time 
at which pinch off would take place for a Newtonian liquid. 
The straight lines are power-law fits yielding an exponent of 
0.66 ± 0.01. The dashed lines are exponential fits to the dis- 
turbance amplitude 1 — h m in/R- The squares in the inset 
are the corresponding growth rates, the line is the theoretical 
prediction for the most unstable Rayleigh mode of an inviscid 
fluid. 



The transition from the symmetric to the asymmet- 
ric shape can be understood by considering the Newto- 
nian motion alone, since the onset of elastic effects is 
very abrupt for low- viscosity solvents. We confirm this 
observation by noting that the motion of the minimum 
neck radius up to just before the onset of elastic effects 
is well described by Newtonian theory, see also 0. For 
the regime before the formation of the filament, at early 
times, the motion is consistent with Rayleigh's theory 
of exponential growth. More quantitatively, this leads 
to the growth of the disturbance of the liquid column 

1/2 

with a rate uir = (0.1187/(/oi? 3 )) , with 7 the surface 
tension and p the liquid density. Corresponding fits to 
hmin are shown in Fig. [3] (dashed lines) , and compared 
to theory in the inset. Fore later times, as the bridge 
becomes more strongly deformed, the well-known simi- 
larity solution 23 h mm = 0.7(j(t ~ i) 2 / p) 1 / 3 for the 
detachment of a low-viscosity droplet takes over. Here to 
is the extrapolated singularity time, at which the droplet 
would pinch off in the absence of polymer. This power 
law corresponds to the straight lines in Fig|3| We see in 
Figl^that the power-law behavior is not observed for the 
smallest nozzle radius: we have not yet transited into the 
similarity solution regime before the elastic forces due to 
the polymer intervene and lead to the formation of the 
filament. 

As the Rayleigh theory is a linear stability argument, 
and the similarity solution is only valid close to the pinch- 
off, there is no theory for the transition between the two 
regimes. However experimental observations show that 
the transition from an asymmetric to a symmetric shape 



occurs at a fixed value of h min /R = 0.17 ± 0.01, corre- 
sponding to a sufficiently strong deformation of the liquid 
bridge. 

Bead formation can now be understood from the in- 
terplay between inviscid dynamics and elastic effects: in 
regions where significant polymer stretching occurs a uni- 
form thread is formed. If elongation rates during the ini- 
tial Rayleigh thinning are sufficiently high for the poly- 
mers to stretch, the symmetric Rayleigh solution trans- 
forms directly in a uniform thread, as shown in Fig. \\]p. 
For larger values of R, however, the elongation rates are 
smaller and polymer stretching sets in later, such that 
the asymmetric similarity solution develops without be- 
ing influenced by the presence of the polymers. 

From Fig^, it is evident that the bead starts to form 
at the bottom of the filament, near the droplet, where 
the fluid neck is thinnest. Since the elongation rate be- 
comes very large in this region, a filament forms in a very 
localized fashion. Meanwhile, the rest of the fluid neck 
continues to evolve as if the fluid were inviscid, eventu- 
ally pinching at the top of the neck. As a result, another 
filament forms, isolating a bead in the middle (cf. Fig. 
^,). Thus the bead is the exact analogue of the so-called 
"satellite drop" 0, which is always formed between the 
main drop and the faucet in the case of a Newtonian 
fluid. However, the polymer prevents pinch-off, and a 
thread remains, connecting the "satellite drop" to both 
the main drop and the faucet. 

If the growth rate ojr of the Rayleigh instability ex- 
ceeds A -1 , such that polymer stretching takes place 
from the start, no bead is formed. Thus R cr u = 
(0.118A 2 7/p) 1//3 « 1mm is the critical radius above 
which a bead is formed, in good agreement with exper- 
iment for low polymer concentrations. The influence of 
polymer concentration, is relatively minor, in agreement 
with this argument: the critical radius changes only by 
a factor of two while the concentration varies by two 
orders of magnitude. The slight increase of R cr it with 
concentration reflects an increase in the relaxation time 
A, perhaps indicating overlap [Tl| . 

The growth of elastic stresses due to the stretching 
of the polymer, can be calculated as follows. The last 
stages of the inviscid motion are very fast, so the poly- 
mer does not have time to relax. In such a case the axial 
polymer stress is given by the total deformation, lead- 
ing to a z — (To/h 4 [12L 1 1 3| . in excellent agreement with 
the results of our numerical simulations. The constant 
(To = r/ p ho/(p\) 13] is determined by the neck radius h Q 
at which the timescale of the flow has become shorter 
than the polymer timescale A, where v v is the polymeric 
contribution to the viscosity. 

For simplicity let us focus on the smallest pipette 
radii, for which we can estimate that ho ~ R; equating 
elastic and surface tension forces, we infer that threads 
start to form when the minimum radius is h t hread/R ~ 
(r]pR/(X"f)) 1 ^ 3 . This means the total deformation the 
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FIG. 4: The migration of the primary bead for the 50ppm, 
R = 1.5mm run (pictures are 1x4 mm). The profiles are 
shown in time steps of 2ms. The bead is pushed upward 
against the direction of gravity. 



polymers undergo before a thread forms is proportional 
to A~ 4 / 3 , where A is proportional to the solvent viscos- 
ity, and thus relatively small for water. Our experiments 
show that the thinning of the filament is exponential, i.e. 
the elongational rate eoj is constant. Exponential thin- 
ning of filaments is generally believed 0, 0, 0] to result 
from a balance of surface tension and elastic forces, con- 
comitant with exponential stretching of polymers. 

The present simulations however do not agree with 
this simple picture: by the time filaments are formed, 
polymer extension is so strong that non-linear effects 
become important within the present FENE-P model- 
ing. We have chosen the finite extensibility parameter 
b = 2.5 x 10 4 , in agreement with [jj and obtain reason- 
able agreement between the observed thinning rates in 
the simulation and the experiments. However, in the sim- 
ulations the thinning rate only remains constant for less 
than a decade, in disagreement with the experimental 
observations. Some experimental confirmation of non- 
linear effects in the constitutive equations is provided for 
by our finding of a weak dependence of eo on the noz- 
zle radius. However, serious questions remain as for the 
applicability of the FENE model to strongly nonlinear 
phenomena in dilute polymer solutions f|15l|). 

After the discussion of secondary bead formation it 
is worth examining the subsequent evolution of beads, 
which can be quite varied. Li and Fontelos |l6| pre- 
dict several scenarios of which we found at least three 
in our experiments: The disappearance of a bead due 
to stretching, bead migration, and fusion of threads due 
to differences in the capillary pressure. Bead stretching 
(between the dashed and full lines in Fig. [2J occurs if 
the bead is small and the pressure difference between the 
bead and the surrounding filament is small. Bead mi- 



gration usually follows gravity, but smaller beads may 
be driven upward by pressure forces (cf. Fig. where 
they fuse with the upper reservoir. Small beads also get 
sucked into larger ones (cf. Fig. ^), because their inter- 
nal pressure is higher. 

In conclusion we have measured and described the de- 
tachment of a drop of a visco-elastic liquid, and deter- 
mined the conditions for the formation of a secondary 
bead in a the detachment process. Our numerical sim- 
ulations describe the drop formation process in very fa- 
vorable agreement with the experiments, and allow us 
explain the qualitative features of the 'phase diagram' 
for the formation of secondary beads. However, for later 
times strong non-linear effects become important, and 
the agreement between experiment and simulation is less 
good. 
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